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Quasi-Quark Spectrum in the Chiral Symmetric Phase 
from the Schwinger-Dyson Equation 

Masayasu Harada,*) Yukio Nemoto**) and Shunji Yoshimoto***) 
Department of Physics, Nagoya University, Nagoya 464-8602, Japan 



We non-perturbatively study the fermion spectrum in the chiral symmetric phase from 
the Schwinger-Dyson equation with the Feynman gauge, in which we perform an analytic 
continuation of the solution on the imaginary time axis to the real time axis with a method 
^ I employing an integral equation. It is shown that the fermion spectrum has two peaks, which 

, correspond to the normal quasi-fermion and the plasmino, although these peaks in the strong 

coupling region are very broad, owing to multiple scatterings with gauge bosons. We find 
' that the thermal mass of the quasi-fermion saturates at some value of the gauge coupling, 

^SJ \ beyond which the thermal (pole) mass satisfies M ^ T, independently of the value of the 

gauge coupling. We also comment on the appearance of a three-peak structure in the fermion 
spectrum as a non-perturbative effect. 



§1. Introduction 



Recent experimental studies carried out at the Relativistic Heavy Ion Collider 
(RHIC) have yielded unexpected findings. Among them, it has been found that 
the collective flow of the created matter behaves like a perfect fluid. This suggests 
that the quark-gluon plasma (QGP) near the phase transition is a strongly inter- 
lO . acting system. This result seems to be consistent with the findings of recent studies 

employing lattice QCD, which show that the lowest charmonium state survives at 
temperatures (T) higher than the critical temperature (Tc).^-* Theoretically, the ex- 
istence of hadronic states in QGP was predicted many years ago on the basis of 
\ symmetry arguments of the chiral phase transition. There are also many recent 

■ studies investigating possible hadronic bound states in QGP which were motivated 

by the RHIC experiments.^-*' 
\ The study of quarks and gluons, as well as hadronic states, is also important, 

^ ■ because they are liberated degrees of freedom in QGP. The rapid change of the energy 

density around Tc obtained in lattice QCD suggests that quarks and gluons actually 
come into play in thermodynamics. However, in general, particles in medium have 
different spectra from those in vacuum, and thus it is quite nontrivial to determine 
whether the quasi-particle picture of quarks and gluons remains valid in strongly 
coupled matter, such as that described by QGP near Tc- It is shown in Ref. 6) that 
the quark spectrum near Tc has two massive modes with thermal (pole) masses of 
the order of T, taking the gluon condensate into account. In Ref. 7), the dispersion 
relations of the quark and the gluon above Tc of the deconfinement transition are 
computed in quenched lattice QCD using the maximum entropy method. The results 
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of those works are that the thermal masses of the quasi-quark and the quasi-gluons 
are larger than T near Tc, while no pure collective modes like the plasmon or the 
plasmino were found. Recently, the thermal mass of the quark was also studied in 
quenched lattice QCD with two-pole fitting for the spectral function. In that work, 
it was found that the thermal mass is of the order of T near Tc- In Ref. 9), with 
a Brueckner-type many-body scheme and data from the heavy quark potential in 
lattice QCD, both the thermal mass and the width of the quasi-quark are found to 
be of the order of T for T = l-2rc. In Ref. 10), it is shown in the Nambu-Jona- 
Lasinio model that the quark spectrum near Tc has three peaks, which result from 
coupling with fluctuations of the chiral condensate, and that the thermal masses of 
the normal quasi-quark and the plasmino are of the order of T. 

In this paper, using a gauge theory, we investigate the fermion spectrum in the 
chiral symmetric phase. In the weak coupling region at finite T, the HTL resumma- 
tion has been established, and the fermion spectrum is well understood at leading 
order. In this study, we compute the fermion spectrum over a wide range of value 
of the coupling constant using the Schwinger-Dyson equation (SDE) with the Feyn- 
man gauge. The SDE used in this paper includes the HTL of the fermion self-energy 
at leading oder, and it thus reproduces the fermion spectrum of HTL in the weak 
coupling region. In the strong coupling region, by contrast, the SDE incorporates 
non-perturbative corrections as an infinite summation of a certain kind of diagrams. 
In addition, one of the advantages of our approach is that the SDE respects the 
chiral symmetry and describes dynamical breaking. This contrasts with the fact 
that it is difficult to respect the chiral symmetry on a lattice. Furthermore, it is 
straightforward to extend our formulation to finite density. 

The contents of this paper are as follows. In ^ we formulate the SDE in the 
imaginary time formalism with the ladder approximation. We perform an analytic 
continuation of the fermion propagator to the real time axis numerically with a 
method employing an integral equation. In f|3l using the retarded fermion propa- 
gator obtained in this way, we investigate the fermion spectrum above Tc, focusing 
on dependences on the gauge coupling. Section U] is devoted to a brief summary 
and discussion. In the appendices, the notation and some technical details of the 
calculations given in the text are presented. 

§2. Schwinger-Dyson equation at finite temperature 

In this section, we introduce the SDE with the ladder approximation at finite 
T. The SDE has been employed in analyses of strong coupling gauge theories in 
vacuum^^) and in the determination of the phase structure of the chiral and color 
superconducting phase transitions at finite T and/or density. ^^-''^^^ 

In the imaginary time formalism, the SDE for the full fermion propagator 
S{iuJn,p) is given by 
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g'^jfj,S{iuJm, k)ryV^"' {iuJn - iuJm,p- k) , (2-1) 



(2vr) 



where F^, and P^'^are the full fermion-gauge boson vertex function and the Matsubara 
propagator for the gauge boson, respectively, and w„ = (2n + l)7rT is the Matsubara 
frequency of the fermion. (For the definition of the Matsubara Green function, see 
Appendix [XI) The quantity g is the fixed gauge coupling. In this study, the current 
fermion mass is taken to be zero, and thus the free fermion propagator is given by 
'^frcc ~ = P ' 1 ~ i^nlo- Here, S is restricted by rotational invariance in space 
and parity invariance to take the form 

with p = \p\- Equation (j2-l|) is not closed for the full fermion propagator. Here, we 
use the ladder approximation in which D^'^ and are replaced with the tree-level 
quantities, ^^^gg and 'j^, respectively. Then, Eq. (|2-ip becomes a closed equation for 

S''^{iuJn,P) - S^^^{iuJn,P) 

= -g^T ^ j--—^-f^S{iiVrn,k)-^^V'^^^^{iujn-ii^rn,p-k), (2-3) 

where the free gauge boson propagator I'^^g is given by 

VZ,{iu^n - iu;m,p- k) = ^ (^g'^-' - (1 - a)^) , (2-4) 

with l^^ = p^^ — k^^ = {iujn — iujm,p — k) and a being the gauge parameter. Here, 
we should mention the gauge fixing of the SDE with the ladder approximation. The 
ladder approximation, expressed by F^, = ju, implies that the vertex renormalization 
factor is unity, i.e. Zi = 1. In QED, the Ward-Takahashi identity implies Zi = Z2, 
where Z2 is the wave function renormalization factor for the fermion. To satisfy 
this equality, Z2 must also be unity. At T = 0, it can be shown that A(pQ,p) = 
C{po,P) = 1) and thus Z2 = 1 is always satisfied if the Landau gauge is adopted. 
For this reason, the Landau gauge fixing is often adopted at T = to study chiral 
dynamics in QCD. At finite T, however, both A and C deviate from unity, owing 
to thermal effects, even in the Landau gauge. Therefore, there is little advantage 
of the Landau gauge at finite T. Here, we use the Feynman gauge in which 
the analytic continuation of the fermion propagator is simple, as discussed in the 
following section. 



*■* In Ref. 17), a non- linear gauge is formulated to satisfy Zi = Z2 in the SDE with the ladder 
approximation at finite T. 
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By making suitable projections, the SDE (j2-3p can be divided into three coupled 
equations for the functions B{iuJn,p), A{iu!n,p), and C{iuJn,p)- After performing the 
three-dimensional angle integral, they are given by 



m=-N 



= 1 + ^ 1 dk 



N ,A 



2txT sr-^ f k A{iuJm,k)E2{iuJn,P',i^JJm,k) 



^^_^^o k'^A'^{iLOm,k) +ujliC'^{iuJm,k) + B'^{iuJm,k) 



(2-6) 



2vrr ^ k'^C{iuJm,k)E3{iujn,p;i(^m,k) 



CiiuOmP) = I + -. y I dk 



m=-N ■ 



Iq k'^A^{iiOm-, k) + Lolfi'^iiu^rri, k) + B'^iiu^rri, k) 

(2-7) 

where the kernels Ei, E2 and E^ are expressed as 

= + + ' (2-8) 



(27r)3 \ 2pk {p + A;)2 + - a;^)2 

(2-9) 

TP r ■ M 5^ ^t^mi (p - + (a;„ - a;m)2 

""^^^""'^'^""''^ = + + • (2-10) 

Because the integrals in Eqs. (|2-5p - (|2-7p are divergent, we introduce the three- 
dimensional ultraviolet cutoff A. We also truncate the infinite sum of the Matsubara 
frequency at a finite number N. However, N can be taken to be sufficiently large so 
that the following results do not depend on it. 

We obtain the functions A{tLi,p), B{lu,p) and C{uj,p) defined at discrete points 
on the imaginary time axis, uj = iuJn, in the complex w-plane. In the present work we 
are interested in the fermion spectrum in the chiral symmetric phase. To specify the 
chiral symmetric phase, we first determine the critical temperature, Tc, of the chiral 
phase transition from these functions and the effective potential. (See Appendix 
[B] for details.) Then, we investigate the fermion spectrum in the chiral symmetric 
phase, in which B = 0. 

Since the fermion spectrum is defined on the real time axis, we need to perform 
an analytic continuation for the solution of Eq. (|2-3p to the fermion propagator on 
the real time axis. Following a method developed in Ref. 15), we perform an analytic 
continuation of the fermion propagator in the imaginary time formalism. An essential 
point of this method is the fact that if the summation of the Matsubara frequency 
in the loop diagram can be done analytically, analytic continuation to the retarded 
function becomes trivial, i.e. iuJn — > Po + ^e- With this goal, we first perform the 
summation of the Matsubara frequency formally, which can be done by expressing 
the fermion and the gauge boson propagators as spectral representations. Then, the 
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analytic continuation can be done easily, and the remaining integrals can be carried 
out numerically. (For details, see Appendix [Cll 

The equation for the retarded fermion propagator obtained in this way is the 
integral equation 



-1 



m=— oo 
d^k 



Sk 

(2^ 



S{iuJm, k) 

PO- Z- iuJr, 



1u(^b{z,p- k) 



+ 9^ 



(27r)2 



dzj^iSnipo - z,k)j„p'^ {z,p- k) 



, Po — z ^ z 

tanh h coth — 

2T 2T 



(2.11) 



where the subscript R indicates the retarded functions, is the spectral function 
for the gauge boson, and S{iuJm,k) is the solution of Eq. (j2-3p . (The definition of 
the retarded Green function and the spectral function are summarized in Appendix 
1X1 ) This method is known to be a reliable method for the analytic continuation in 
condensed matter physics and is suitable for the SDE, because Eq. (|2-3p is also an 
integral equation that is similar to Eq. (j2-lip . 

It is noted here that Eq. (j2-lip is valid for any gauge fixing. The gauge depen- 
dence is in the form of the spectral function, p'^ , only. In this paper, we adopt the 
Feynman gauge, in which p^ is expressed by 

= -g^-e{po)5{pl-p') , (2-12) 

where e(po) is the sign function, e(po) = il foi" Po ^ 0. For the Landau gauge, on 
the other hand, the evaluation of the double-pole terms in the propagator is rather 
tedious. This is one reason that we adopt the Feynman gauge in this study. 
In the chiral symmetric phase, the retarded fermion propagator is written 

-iSr{pq,p) = ^ + ^777 ^> (2T3) 

D+{po,P) D^{po,p) 



where p = p/p, and D± {pQ , p) is defined as 



D±{po,p) =PoCr{pq,p) ^pAr{pq,p) . (2-14) 

Here, Ar{pq,p) and Cr{pq,p) are the retarded functions obtained from A{iujniP) 
and C{iuJn,p) by analytic continuation, respectively. 

Performing the angular integrals in Eq. (j2Tip and the projection for the Dirac 
index, we obtain the following self-consistent coupled equations for Ar{pq,p) and 
Cr{po,p)- 
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2ttT ^ k'^A{iuJm,p)E2{po,p;iu;m,k) 



Ar{po,p) = 1 + ^ dk 

^ m=-oo •^'^ 



g2 roo r\p+ko-po\ 

+ ic 9 ^ / dko / dk 
167rV J-oo J\p-ko+po\ 

^ k{p^ + fc2 - (po - kof)Aiiipo,p)X{po, ko;T) 

k^Cl{ko,k)-k^Al{ko,k) 

2ttT ^ /"^ k'^C{iuJm,p)E3{po,p;iuJm,k) 



2ttT i 

Cr[po,p) = i + — V / dk 



k'^A^{iuJm-, k) + ujlfi'^iiu^rn, k) 



(2-15) 



+ r dko t"^"'""'"^ dk kkoCR{po,p)X{po,ko;T) 

SvrVo i-oo ° J\p-ko+po\ k^Cliko, k) - k^Aj^{ko, k) ' 

where X(po, ^o; T) is expressed as 

X(po, ko;T)= tanh ^ + coth . (2-17) 

In the above expressions, we have replaced the integral variable zhy k^ = pq — z. We 
should note that A{iuJn,p) and C{iu!n,p) on the right-hand sides of Eqs. ()2-15p and 
(j2-16p are the solutions of the SDE (j2-3p obtained in the imaginary time formalism. 



§3. Spectral function 



3.1. Coupling dependence 

In this subsection, we study the fermion spectrum over a wide range of value 
of the gauge coupling, focusing on the peak structures. Because the temperature 
T is the only infrared scale in the chiral symmetric phase, dimensional quantities 
below are scaled by T. In the following, we fix the cutoff A as T/A = 0.3. We have 
confirmed that the present results depend very little on the choice of A. (See the 
end of this section.) 

The fermion spectral function is given by 

P±{Po,p) = --Im— — , 

TT D±{po,p) 

where (p-) denotes the fermion (anti- fermion) sector. The form of (/0_) can 
have several peaks, reflecting collective excitations, e.g. the normal quasi-fermion 
(quasi-antifermion) and the anti-plasmino (plasmino) in the HTL approximation. 
We consider only in the following, unless otherwise stated, noting the relation 

P-{Po,p) = P+{-Po,p)- 

We first investigate the fermion spectrum in the weak coupling region. The 
spectral function for p = and g = 0.115 is shown in Fig. [TJ We see that 
p+ has two sharp peaks, one in the positive energy region, corresponding to the 
normal quasi-fermion, and the other in the negative energy region, corresponding to 
the anti-plasmino. The peak positions of /?+ and p_ for positive energy and finite 



(3-1) 
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Fig. 1. The spectral function p+ at rest (p — 0) for g — 0.115. The sharp peak in thepo/T > (< 0) 
region corresponds to the normal quasi- fermion (the anti-plasmino) . 
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Fig. 2. Peak positions of the spectral functions for the normal quasi-fermion (circles) and the plas- 
mino (squares) for g = 0.115. These are determined from the peak positions of the spectral func- 
tion. The solid and the dashed curves represent the dispersion relations for the quasi-fermion 
and the plasmino in the HTL approximation, respectively. The thin solid line corresponding to 
po ~ p is also shown for convenience. 



momentum are plotted in Fig. [21 which displays approximate dispersion relations of 
these excitations. We see that the dispersion relations obtained from the SDE are 
in good agreement with those obtained with the HTL approximation, which is valid 
in the weak coupling region. The narrow widths of the peaks imply that the spectra 
from the SDE are quite similar to those in the HTL approximation in the weak 
coupling region. This result is natural, because the SDE here contains the fermion 
self-energy in the HTL approximation, which is dominant in the weak coupling region 
in our approach, as expected. 

For the case of strong coupling, we plot the spectral functions p-^ for g = 1.15,3 
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Fig. 3. The spectral functions p+ for g — 1.15 (left), g = 3 (middle), and g — 6.9 (right). The 
peak for po > (< 0) corresponds to the normal quasi-fermion (the anti-plasmino) . The figure 
is clipped at p+T = 0.3. 
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Fig. 4. The spectral functions p+ at rest (p = 0) for g = ?>. The solid curve represents the spectral 
function obtained from the SDE, the dashed curve that from the one-loop self-energy, and the 
dotted line that from the HTL approximation. 

and 6.9 in Fig. [3l All the spectral functions appearing in Fig. [3] have two peaks 
but with broad widths. This shows that the normal quasi-fermion and the plasmino 
appear even in the strong coupling region. We see that as the coupling becomes 
larger, the widths become broader. In the case of extremely strong coupling, e.g. 
g = 6.9, the peaks are so broad at low momenta that the quasi-particle picture is no 
longer valid. For these three couplings, we have confirmed that as the momentum 
becomes higher, the peak of the plasmino becomes smaller, and the peak of the nor- 
mal quasi-fermion becomes larger and sharper. This behavior is due to the fact that 
thermal effects become smaller as the momentum increases, and thus the spectrum 
approaches that for the free fermion. This momentum dependence of the heights of 
these peaks is similar to that of the residues of the peaks in the HTL approximation. 

In the left and middle panels of Fig. [Sj there is another small peak at the origin 
in addition to the two obvious peaks. We comment on this peak in i ^3.2[ 

In Fig. m we plot the spectral function p+ for 5 = 3 and p = in the positive 
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Fig. 5. Coupling dependence of the thermal mass. The solid curve represents the thermal mass 
obtained from the SDE, and the dashed line that from the HTL approximation. 




Fig. 6. Cutoff dependence of the ratio M/T for g = 6.9. The solid line represents M/T = 1.19. 



energy region and compare it with that obtained from the one-loop self-energ}El] and 
that obtained in the HTL approximation. This figure shows that the peak obtained 
at one-loop order is already broad in the strong coupling region. It also shows that 
the peak from the SDE becomes even broader owing to non-perturbative effects. 
Because the SDE contains effects of multiple scatterings with gauge bosons through 
the self-consistency condition, the broadening of the peaks could be understood from 
the fact that the probability of gauge boson emission and absorption from a fermion 
increases rapidly with the coupling. 

Now, we interpret the peak positions for p = as the thermal masses of the 



*■* Here, the one-loop self-energy includes the T-independent part evaluated with the cutoff A. 
In Ref. 19), on the basis of the one- loop self-energy without the T-independent part, it is shown 
that the thermal mass is slightly larger than that in the HTL approximation, with the difference 
being less than 5%. 
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Fig. 7. The spectral functions p+ at rest [p = 0) for g = 3 (dashed curve) and g = 2 (solid curve). 

excitations irrespectively of the widths. In Fig. [5] we plot the coupling dependence of 
the thermal mass together with that obtained in the HTL approximation. We note 
that the thermal masses of the normal quasi-fermion and the plasmino are always 
the same in the chiral limit at zero density. ^'^^ We see that in the weak coupling 
region, the thermal mass from the SDE almost coincides with that in the HTL 
approximation, and both of them are proportional to the coupling, g. In the strong 
coupling region, by contrast, the thermal mass from the SDE becomes saturated at 
M ~ T and is almost independent of the g. 

Before ending this subsection, we study the dependence of the resultant fermion 
spectrum on the choice of the ultraviolet cutoff A. In Fig. [6] we show the A- 
dependence of M/T for g = 6.9. We see that the ratio is almost constant, M/T ~ 1.2, 
for A/T > 2.5. This implies that the thermal mass is determined by only T which 
gives an infrared scale of the system. Furthermore, we have confirmed that the 
fermion spectral functions shown in Fig. [3] have little dependence on the cutoff for 
Pq ~ M. The shapes of the spectral function in the large pq region are slightly 
modified when we use smaller cutoff, while we see the tendency for the peak of the 
normal quasi-fermion becomes larger and sharper as the momentum becomes higher. 
Thus, we conclude that the thermal mass M, as well as the tendency of the shape, 
is not affected by cutoff artifacts for A/T > 2.5. For A/T < 2.5, however, deviation 
from a constant value is seen. There cutoff effects might start to contribute. 

3.2. Comment on the three-peak structure 

Finally, we give a comment on the peak around the origin that appears in the 
spectral function for g = 3, as seen in the previous section. This peak is more 
pronounced for weaker coupling, = 2, as shown in Fig. [TI*^ For stronger coupling, 
such a peak is not seen. One possible reason for this behavior is that large widths 

*' In Fig. [T] the spectral function near po/T = is not plotted because of the infrared cutoff 
introduced in the numerical analysis. A sharp peak with a width narrower than the infrared cutoff 
should exist. In the present analysis, we cannot confirm its existence due to numerical error. 
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of all the peaks due to multiple scatterings with gauge bosons may conceal the peak 
structure around the origin. 

It has been shown in Yukawa models a similar peak appears in the fermion 
spectral function owing to the coupling with a massive boson at finite T.^^)'^^) In 
this case, the mass of the boson is essential for the formation of the three-peak 
structure in the fermion spectral function. In our case, however, the fermion couples 
to a massless gauge boson, and thus another mechanism makes it possible to form 
the three-peak structure. Because the SDE is a self-consistent equation for the 
full fermion propagator, the internal fermion in the self-energy has a thermal mass. 
Therefore, the system is similar to that consisting of a fermion with a thermal mass 
coupled to a massless gauge boson. Here, as an example of such a system, we consider 
an interaction of a fermion with a thermal mass obtained in the HTL approximation 
and a massless gauge boson at one loop. We summarize the details of this analysis 
in Appendix [Dj Figures [10] and [11] in Appendix [D] plot the imaginary part of the 
fermion self-energy and the spectral function at rest in this model. We see that the 
imaginary part of the self-energy possesses two peaks, which lead to the three-peak 
structure in the spectral function, as shown in Refs. 21) and 22). Therefore, the 
existence of the three-peak structure of the spectral function from the SDE can be 
understood from this result. 

§4. Summary and discussion 

We have investigated the fermion spectrum in the chiral symmetric phase em- 
ploying the Schwinger-Dyson equation (SDE) for the fermion propagator with the 
ladder approximation and fixed gauge coupling, in which the free gauge boson prop- 
agator in the Feynman gauge is used. Because we employed the imaginary time for- 
malism to construct the SDE, analytic continuation to the real time axis is needed 
to obtain the fermion spectrum. This was done by solving an integral equation for 
the retarded fermion propagator, following a method developed in Ref. 15). This 
method is known as a reliable way to carry out numerical analytic continuations 
in condensed matter physics, and we have confirmed in this study that it is also 
effective even in a relativistic system. 

We have investigated the fermion spectrum in the chiral symmetric phase over a 
wide range of values of the gauge coupling using the method of analytic continuation 
mentioned above. We found that in the weak coupling region (g <C 1), the fermion 
spectrum is quite similar to that obtained in the hard thermal loop (HTL) approxi- 
mation: There exist a normal quasi-fermion and a plasmino whose thermal masses are 
almost the same as those found in the HTL approximation, M ~ Mhtl = gT/ ^/8, 
and their widths are very narrow. (These widths vanish in the HTL approximation.) 
This result is natural, because the SDE here contains the HTL of the fermion self- 
energy, and thus the HTL is dominant in the weak coupling region in our approach, 
as expected. 

In the strong coupling region (g > 1), we find that there exist a normal quasi- 
fermion and a plasmino, as in the weak coupling region. However, their spectra 
deviate from those in the HTL approximation: Both the excitation modes have 
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smaller thermal masses and much broader widths in the low-momentum region. The 
thermal mass M saturates at some value of the gauge coupling, and satisfies M r^T 
independently of the value of the coupling. Although we do not have any plausible 
explanation for the saturation of the thermal mass at present, the saturated value 
M ~ T is consistent with the results of other model analyses.^^'^^"^'^)'^^) In Ref. 9), 
the broad width is described in terms of the imaginary part of the self-energy, while 
in our study, it is described in terms of the spectral function. Furthermore, we 
have shown that as the gauge coupling becomes larger, the widths rapidly become 
broader. This indicates that the broadening of the widths can be understood from 
the fact that the probability of gauge boson emission and absorption from a fermion 
increases rapidly with the coupling. 

Our result shows that as the momentum p becomes higher, the peak of the 
plasmino becomes smaller and the peak of the normal quasi-fcrmion becomes higher 
and sharper. Thus the spectrum approaches the free-fermion one. This behavior is 
similar to that in the weak coupling region and is natural because thermal effects 
are weaker in general for p > T. The free-fermion spectrum in the high momentum 
region is consistent with the quark recombination or coalescence picture,^^) in which 
the quasi-quark picture is valid. 

One reason we have adopted the Feynman gauge is for the simplicity of the 
numerical analytic continuation mentioned above. We expect that the qualitative 
features of the fermion spectrum discussed in this paper do not change even if we 
adopt other gauge fixings. 

We now discuss what our result can provide regarding the quark spectrum in 
the real-life QCD at high temperature. In the very high temperature region of 
QCD, the HTL approximation is valid for studying the quark spectrum. In the HTL 
approximation, the fixed coupling, which is the running coupling at an energy scale 
of order T, say g{T), is used, because the quarks and gluons with energy T give the 
dominant contribution. In studying the quark spectrum in the lower T region, the 
diagrams other than the HTL start to give non-negligible corrections, since the (fixed) 
coupling g{T) is larger. For this reason, we can state that, in our approach, we include 
a part of the corrections from all the ladder diagrams by solving the Schwinger-Dyson 
equation. Actually, our result reproduces that of the HTL approximation in the weak 
coupling region, i.e. in the very high temperature region where g{T) is small. In 
the lower temperature region, where g{T) is large, our result in the strong coupling 
region suggests that the corrections from the ladder diagrams have the following 
effects on the qualitative structure of the quark spectrum: (1) The thermal masses 
of quark and plasmino are not proportional to the gauge coupUng g{T), becoming 
M T, and (2) their widths are very broad, due to multiple scatterings with gluons. 

In this study, we have used the free gauge boson propagator in the SDE. In 
medium, however, the gauge boson spectrum may also change significantly, as in the 
HTL approximation. One possible extension to include such an effect is to use the 
gauge boson propagator with Debye screening and a dynamical screening that has 
the same form as that in the HTL approximation. Such a study is left as a future 
project. 
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Appendix A 

Notation 



Here, we briefly give definitions of some of the quantities used in this paper. 
The Matsubara Green function is defined by 

S{iun,p) = F.T.(r^V(r,fOV^(r',P))T , (A-1) 

where (• • • )t stands for the thermal average with temperature T. The retarded 
Green function for a fermion Sr and the spectral function for a fermion are defined 

by 

SRit,f;t',r') ^ e{t-t'){{i^it,f)Mt',f')})T , (A-2) 
piPo,p) = ^F.T.{{i;{t,f),lp{t',r')})T • (A-3) 

The retarded Green function and the Matsubara Green function are related to the 
spectral function by the spectral representations as 

S{zcOn,p) = - r dp',P^^^ , (A.4) 

iSRipo,p} = - r dp', . (A.5) 

These functions for gauge bosons can be written in the same forms. 

Appendix B 

Phase Structure 



In this appendix, we compute the critical temperature Tc to specify the chiral 
symmetric phase. We introduce the effective potential for the fermion propagator to 
determine the phase structure and Tc,^^^ 

ViS) = T / 7^{lndet[5(za;„,^]+tr[^5(za;„,^]} 

n=— oo ■ ^ ' 

^ 1„27.2 f d^P f d^k 

^ r ^ J (27r)3 J 



(27r)3 

n,m=—oo " / " V / 



X tr ^S{iuJn,p)'yij.S{iuJm,k)'yu 'D^"'{iujn - iuim,P- k) , (B-1) 
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Fig. 8. Temperature dependence of the function B{iun = Q,p = 0)//! (left) and the effective 
potential Kci(r)/yl'* (right) for g = 6.9. 



where In, det and tr are taken in the spinor spaces. The SDE is obtained as the 
stationary condition for the effective potential ()B-ip . Note that the value of the po- 
tential is meaningful only at the stationary point; i.e. the potential can be evaluated 
only when the solution of the SDE is substituted. There always exists a solution of 
the SDE with 5 = 0. In the low temperature region, we also have a solution with 
B ^ 0, which corresponds to dynamical chiral symmetry breaking. Let Sw denote 
the solution with B = and Sng the solution with B ^ 0. They are parameterized 
as 

[Sw{ii^n,p)r^ = Aw{iu;n,p)p- 7 - Ci4/(iu;„,p)iu;„7o , (B-2) 

To determine which of the solutions Sw and Sng gives the real vacuum, we compare 
the values of the effective potential with each solution by computing the following 
difference: 



Kei(r) = viSNG) - vis 



W 



^ oo „A I 

-^T dppHln 



p^A%Q{iuJn,p) +ujiC$fQ{iuJn,p) +B%Q{iu:n,p)' 



P^ Al^{iLVn,p) + UjlC^iiLVn^p) 
^ p^ANG{i^n,p) +(^nCNG{i^n,p) 



P^A%Q{iuJn,p) +ujlC'^Q{iuJn,p) + Bj^Q{iuJn,p) 

p'^Aw{iuJn,p) + ujlCw{iuJn,p) \ 
P^Al^{iuJn,p) + ujlC'^{iuJn,p) I ' 



(B-4) 



For V^e\{T) < 0, the chiral broken vacuum is realized, while the chiral symmetric 
phase becomes the true vacuum for V^eiiT) > 0. We determine Tc as the smallest 
value of T at which Vj-ai{T) vanishes. 

We plot the T dependence of the function B{iuJn,p) for the smallest Matsubara 
frequency at p = in the left panel and that of the potential V^ei in the right panel of 
Fig. [HI It is seen that the function B and the potential T^ei simultaneously go to zero 
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g 

Fig. 9. Coupling dependence of the critical temperature Tc/ A. The chiral symmetric (broken) 
phase is above (below) the line. 

smoothly with T. This imphes that the chiral phase transition is of second order. 
The value of the critical temperature is determined as Tc/A ~ 0.18 for g = 6.9. The 
value Tc/A depends on the coupling. In Fig. [9l we plot Tc/A for each coupling. Note 
that the chiral broken vacuum is not realized even at T = for g < gc — 5.1. 



Appendix C 

Analytic Continuation for Solutions 

of the Schwinger-Dyson Equation 



In this appendix, we explain the details of the method for the analytic continu- 
ation used in the present analysis, following Ref. 15). 

Let us start with the spectral representations for the fermion and gauge boson 
propagators. 



oo 



V^y{iujn - iu}m,p- k) = - I dz-—^ — -. , (C-2) 



oo 



where pp and represent the spectral functions of fermion and the gauge boson, 
respectively. The function Sii{z, k) is the retarded fermion propagator, which is 
analytic in the upper half complex plane of the variable z. 

Using the above spectral representations, the summation over the Matsubara 
frequencies in the SDE (|2-3p is performed by replacing the sums over frequencies 
with the following contour integral: 

"S^ fiko = iuJrn) = <f ^/(^o)^tanh^ . (C-3) 



m=— cxD 



Jc 27ri ^ ^2 2T 
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Then we replace the Matsubara frequency as icOn Po + i^- Thus we obtain the 
SDE for the retarded fermion propagator as follows: 



'R \FU,y; "^ijfrce' 
g2 f (fk 



vr J (27r)3 

lm[iSR{z,k)] 



dz 



dz' 



7/. 



Pq — z — z' + ie 



K - tanh — - + coth — - 
^ 2 V 2r 2T ) 



(C-4) 



By using the relations Re[/(z)] = Im[i/(z)] and Im[/(2:)] = Re[— ^/(z)], Eq. ()C-4p 
can be rewritten as 

iS-^{pQ,p) - iS^^^^J^PQ^p) 



2 roo 



/ dz' 



4tt 
X < Im 



d^k 



dziSR{z, k) 

( . ^ . z'\ 

tanh—- + coth— - 

V 2r 2T J 

dziSji{z, k) 



-, + 

Po — z — z' + le po — z — z' — le 



iRe 



Po — z — z' + ie Po — z — z' — ie 



tanh—- + coth—- >7jy 
V 2T 2T J y 



(C-5) 



Because the retarded fermion propagator, Sr, should be analytic in the upper half 
complex z plane, the integrand of the right-hand side can have poles only at z = 
Po — z' + ie and z = i{2n + l)7rT for n > 0. Then, the integration over z can be 
replaced by the contour integral. Then, performing the contour integral over z, we 
obtain 

iSR\po,P)-iS]i\eeiPO^'P) 



27r 



d^k 



X < Im 



(27r)3 

oo 

2iTTT ^ iSR{iuJm,k) 



+ 



m=0 

iTTiSR{po — z' , k) ( tanh 



Po — iojm — z' + ie Po — iuJm — z' — ie 



zRe 



2i7rr ^ iSR{iujm,k) 



m=0 



^'^-T^; 1- coth— - 

2T 2T 



1 



1 



Po - iuJm - z' + ie Po- ioJm - z' -ie^ 
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i-KiSiiipo - z', k) ftanh ^"^^^ + coth^ 



(C-6) 



In the above expression, we can drop the term ie, because the Matsubara frequency, 
LVmi is always non-zero. Using the relations lva.f{z)—iR.ef{z) = —if{z) and lm.f{z) + 
iRef{z) = if{z)*, we obtain 



iSji^{po,p) - iSii\^^{po,p) 
(fk 



dz' 



p%{z ,p- k)jf. 



S{iuJm,k) S*{iu;m,k) 



+ 



m=0 



PO - ILOm - Z' Po- lUJn 



iSnipo -z'^k)^ ^tanh ^°^^^ + coth^ 



(C-7) 



Finally, because the Matsubara Green function is real, i.e. S{iujrn,k) = S{iujrn,ky 
we arrive at 



iSji^iPo,P) - iSji)^^^{po,p) 
(fik 




dzT ^ 7^ 



m=—QO 



S{iLVm,k) 
Po- Z- ilVr, 



1 



dz-ff,iSR{po - z, k)j^p'^{z,p- k) 



tanh^^^^— h coth— - 

2T 2T 



(C-8) 



With the solution of Eq. ()2-3p substituted into S{iuJm,k), Eq. ()C-8P becomes a self- 
consistent equation for the retarded fermion propagator. In this way, we can perform 
the analytic continuation for the solution of the SDE by solving the integral equation 

del). 

Note that the first term in Eq. ()C-8P is obtained through the naive replacement 
iuJn ^ Po + ic in Eq. (j2-3p . which is easily checked by integrating over z in Eq. (|C-8p . 
The second term appears from the fact that the Matsubara summation must be done 
before the analytic continuation. We can also verify that Eq. (IC-8P returns to the 
original SDE ()2-3p through the replacement po — > iun- 



Appendix D 

Three-Peak Structure in a Simple Model 

In this appendix, we show that the three-peak structure of the fermion spectrum 
obtained in the SDE can be understood from the interaction of a fermion with a 
thermal mass and a free gauge boson. For simplicity, let us consider a one-loop 
diagram consisting of a fermion propagator obtained in the HTL approximation and 
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a massless gauge boson propagator. Then, the fermion self-energy is given by 

U{iujn,p) = -g^T^^ J -^^-i^,S^'^^{ii^n,,^-iuV'^^^^{iuJn-iuJm,p-k) , (D-1) 

where 5^^^ is the fermion propagator obtained in the HTL approximation and X^f^gg 
is the free gauge boson propagator in the Feynman gauge. The retarded fermion 
propagator in the HTL approximation is decomposed as 



(D.2) 



where is for the fermion and for the anti-fermion. The quantities 

are related to the spectral function as in Eq. (j3-ip . The spectral functions 
can be written as 



pI^^{jPo,p) = Z±{p)6{po - oo^ip)) + Z^{p)6{po + uj^{p)) + pUpo,p) , (D-3) 



where Z±{k) = {u!^{k) — k^)/2M^r^-^ are pole residues, with M^rp^ = g^T^ /S, and 
are the continuum parts, given by 

/'±(po,p) = ^m2tl(1tx)0(1-x2) 



ipil ^x)± [(I T In 

2p 



x + l 



X - I 



-1 



(D.4) 



with X = po/p. In the following, we limit ourselves to the rest frame = 0). Using 
p^^, we obtain 



Z'o(iu;n,0) = -tr[7oZ'(fcj„,0)] 



p'i''^{z,k)+p^'^{z,k) f^^, pB{z',k) 



(D-5) 



with pb{po,p) = ^{Po)^{Po —p^)- By carrying out the Matsubara summation and 
the replacement iujn — > Po + i^-, becomes 



^o(Po,0) = g^ 



d^k 
(2^ 



dz I dz 



, 1- f{z)+n{z') 
Po — z — z' + ie 



X [pl^'^{z,k)+p'l^'^iz,k)] pBiz',k) 



(D.6) 



with f{E) = [exp(£'/T) + 1] ^ and n{E) = [exp(£'/T) — 1] ^. The imaginary part 
of Uo can be easily calculated as 



Imi^o = |- 
47r 



^^Af{po-k)-n{k)-l) 
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kZ^{k) 



+ 



+ 



kZ^{k) 



kZ+{k) 



kZ+{k) 



+ H '\/, (/(P0-fc)-n(fc)-l) 
(f{po-k)-nik)-l) 
U{po-k)-n{k)-l) 
Xf{Po + k) + n{k)) 
Af{po + k) + n{k)) 
Af{po + k) + n{k)) 

Xf{po + k) + n{k)) 



kZ^{k) 



1 + UJ'_ 

kZ_{k) 



1-UJ'_ 

kZ+{k) 



Po—k+LU-{k)=0 

P(,~k+u)j^{k)=0 
PO+k—uj+(k)=0 
PO+k+ui^{k)=0 
po+k—ui- {k)=0 



2 



+ — I dkk 

Ait 



po+k+LO+{k)=0 

ifiPo -k)- n{k) - l)(p^+(po -k,k)+ pUpo -k,k)) 
-{f{p^ + k) + n{k)){pX{p^ + k,k)+p^_{po + k,k))\ , (D-7) 



with uj'j^ = duj±{k) / dk. 

It is noted that Eq. ()D-7p is divergence-free. The real part of U^, by contrast, 
has an ultraviolet divergence, and here we regularize it using the relation 



Rei7o(po,0) 



vr 



Po 



P'o 



(D. 



with the three-momentum cutoff A. We choose A to satisfy T/A = 0.3. Here, P 
represents the principal- value operator. By using Sq, the spectral function at rest is 
obtained as 



P±(Po,0) 



— Im — 

vr Po 



1 



(D-9) 



^o{Po,0) 

Now, we present our numerical result for g = 2. In Fig. [10] we plot the po de- 
pendence of the imaginary part of the self-energy, ImZ'o(pO)O). We see two distinct 
peaks for |po/-^htl| < Ij which come from the Landau damping and the kinemat- 
ics: The Landau damping implies a non-vanishing imaginary part that is even for 
|po/-^htl| < 1) while the imaginary part vanishes at the origin, because of the sup- 
pression of the decay rate for distribution functions. These two peaks in Imi7o(pO) 0) 
are similar to those in the system studied in Ref. 21) with a massless fermion and 
a massive boson that result from a Yukawa coupling. Because the current fermion 
mass is taken to be zero, the fermion spectrum is determined solely from the struc- 
ture of the imaginary part of the self-energy, as shown in Eq. (|D-9p . Therefore, the 
fermion spectrum should also be similar to that obtained in Refs. 21) and 22), and 
actually form a three-peak structure, as shown in Fig. \TT\ which is qualitatively the 
same as that in Refs. 21) and 22). 
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o 




Fig. 10. The imaginary part of the self-energy for g — 2, T/A = 0.3 at rest (p = 0). 



+ 

Q. 




Fig. 11. The spectral function for g = 2, T/A = 0.3 at rest (p = 0). 

The peak around the origin in Fig. [11] is very sharp, because the imaginary 
part of the self-energy vanishes at the origin, as mentioned above. Contrastingly, 
the peak at the origin obtained using the SDE is not so sharp, due to the non- 
perturbative iteration: If we use the spectral function ()D-9P for the internal fermion 
line in Eq. ()D-ip instead of 5htl as an effect of the iteration, the imaginary part of 
the resultant self-energy will be non-zero even at the origin, and the peak around 
the origin should become broader. 
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